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Over time, there has been a continuous surge in the demand for electrical 
energy, necessitating the development of larger and more intricate electrical 
power networks. These extensive networks pose a significant challenge, 
primarily in the form of considerable loss of electrical energy, which, if not 
effectively addressed, may lead to persistent and imperceptible losses. In 


response to this challenge, this research proposes the application of quantum 
binary particle swarm optimization (QBPSO) for the coordinated 
management of on-load tap changers (OLTC) in loaded transformers within 
a distribution network, with a specific emphasis on reducing power losses. 
The experimental results demonstrate that the implementation of QBPSO 
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1. INTRODUCTION 

With the passage of time, the demand for electrical energy continues to surge. This escalating need 
necessitates larger and more complex electrical power networks. However, as these extensive networks 
evolve, they bring forth a significant challenge. One such challenge is the substantial loss of electrical energy 
within the network, which, if left unaddressed, can lead to ongoing and imperceptible losses. 

The growing demand for electricity due to population expansion has necessitated the installation of 
power systems with a larger capacity. However, as the capacity of a power system increases, so does its 
power loss and operational cost [1]-[3]. Numerous approaches have been proposed to reduce power loss, 
such as preventing voltage drop, reducing unbalanced phase loads, correcting power factors using capacitor 
banks, and optimizing the loading of power transformers [4]-[6]. Photovoltaics (PVs) have also been 
integrated into power systems to control voltage level fluctuations and reduce power loss [7]-[9]. Static 
compensator distribution [10] and voltage control on the transformer using on-load tap changers (OLTC) can 
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also be used to achieve the same objective [11]. Distributed generators (DGs), wind turbines, and PVs can 
additionally support extra power demand and stabilize supplied voltages [12]-[16]. However, OLTC 
coordination and the exchange of reactive power from DGs, wind turbines, and PVs are time-consuming and 
expensive to install and maintain. 

In this research, we propose the utilization of quantum binary particle swarm optimization (QBPSO) 
to facilitate OLTC coordination and minimize power losses within a distribution network. The choice of 
particle swarm optimization (PSO) as the optimization technique is rooted in its advantageous qualities, 
including speed and efficiency, ease of implementation, scalability, adaptiveness, and a balanced exploration- 
exploitation trade-off. This paper delves into issues related to power losses within distribution systems and 
elaborates on the proposed approach for OLTC coordination. Experimental results stemming from the 
implementation of the suggested method are presented, and a concise conclusion is provided. 


2. METHOD 
2.1. Transformer 

The electromagnetic induction principle is utilized by an electrical device known as a transformer to 
transfer alternating voltage from one level to others [17]. The transformer comprises two types of coils, the 
primary and secondary coils, which are electrically separated but magnetically connected, as depicted in 
Figure 1. The working principle of a transformer is based on Ampere’s law, which states that an electrical 
current can generate a magnetic field, and Faraday’s law, which states that a magnetic field can generate an 
electric current [18], [19]. When the primary coil is connected to an alternating voltage supply, alternating 
flux is generated inside the laminated core. This flux induces self-induction in the primary coil and mutual 
induction in the secondary coil, resulting in magnetic flux in the secondary coil. If the load is connected to 
the secondary coil, the secondary current flows. 

A simple transformer is composed of two insulated coils or wires and an iron core. The ratio of turns 
on a transformer can be determined using (1) and the number of turns on the transformer can impact the 
voltage and current produced on the secondary side. This concept is related to the calculation of the number 
of transformers turns using (2), where Np and Ns represent the number of primary and secondary turns, and 
Vp and Vs are the voltages on the primary and secondary sides, respectively. Lastly, Ip and Is indicate the 
currents on the primary and secondary sides. 


Ns 
aa (1) 
Ns V, Is 
= Bats (2) 


Figure 1. Transformer circuit 


2.2. Problem formulation in the radial distribution network 

The primary objective of this paper is to enhance the efficiency and economic viability of a 
distribution system by minimizing power loss while ensuring that the output voltage remains within specified 
limits. The overarching goal is to optimize the overall performance of the system by addressing both real and 
reactive power losses. Real power, which represents the actual energy consumed or dissipated as heat, and 
reactive power, associated with the phase difference between voltage and current in AC circuits, are targeted 
for reduction. In (3) and (4), as referenced from sources [20], [21], play a pivotal role in guiding the 
minimization process. By focusing on these equations, the paper aims to provide a systematic approach to 
achieving optimal power management, thereby contributing to the economic efficiency and sustainability of 
the distribution system. 
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2.3. Power flow in the distribution system 

The power flow solution in a distribution system diverges from that in a transmission system due to 
the radial-connected network characteristic of distribution systems. A method commonly employed for 
power flow calculations in radial distribution systems is the bus injection to branch current—branch current to 
bus voltage (BIBC-BCBV) method, as outlined in references [22], [23]. This technique facilitates the 
accurate determination of power flow dynamics in the context of radial distribution networks. Figure 2 allows 
for the arrangement of the bus injection to bus current (BIBC) matrix, which can then be used to simplify (5) 
into (6). 


[B, B> Bz B, Bs] =[1111101111001100001000001 Jl, I iy Is Ie] (5) 


[B] = [BIBC][/] (6) 


Figure 2. Single line diagram radial distribution network 


The branch current to bus voltage (BCBV) matrix can be formulated from the voltage equation, as shown 
in (7). 


[Vi — Vz V, — V3 Vı — V4 Vı — V5 V1 — Ve ] = 

[Z12 0 0 0 0 Z12 Z23 0 0 0 Z12 Z23 Z34 0 0 Z12 Z23 Z34 Zas 0 Z12 Z23 0 0 Z36 | 

[B, B, B3 B, Bs | (7) 
Similarly, (7) can be simplified to (8). 

[AV] = [BCBV][B] (8) 
If (6) is substituted into (8), then AV can be expressed as (9) which in turn can be simplified to (10): 

[AV ]=[BCB V][BIBC][T] (9) 

[AV ]=[DLF]I1] (10) 


The complete power change can be obtained by iterating (5) through (13). 


i(k) = (E (11) 
[AVk] = [DLF][Ik] (12) 
[Vk +1] = [V1] - [AVK] (13) 
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2.4. Quantum binary particle swarm optimization 

The PSO algorithm, developed by Kennedy and Eberhart in 1995, draws inspiration from bird and 
fish colony intelligence [22]-[25]. In (14) and (15) are representative of the standard PSO algorithm. The 
introduction of the inertia weight as a diversity controller in the original PSO modifies the particle update 
equation. 


vi(t +1) = vit) +c. rand (pp —x,(t)) + cz. rand (Pg — x,(t)) 14) 
xit +1) = x(t) +0,(t +1) (15) 


The standard procedure for implementing the PSO algorithm includes initializing the velocity, 
position, and PSO parameters, updating the velocity and position of the particles using (14) and (15), 
evaluating the fitness function with Pp and Pg update, comparing each candidate Pg from the fitness 
function value to get the best Pg value, and returning to step 2 if the value of Pg is not the best value. The 
iteration limit is performed to get Pg with the least or the best fitness function [26], [27]. The binary particle 
swarm optimization (BPSO) algorithm updates Pi and Pg in the swarm, with a difference in the 
interpretation of the velocity compared to the PSO standard. The range of the velocity in BPSO is limited to 
[0,1]. In (16) is used to obtain the new position of the particle, while its velocity is determined by (17). 


v'i (t) = sig((t)) = ay n 
$ -{tifry < sig(v(t + 1)) 
Xij (t + 1) { 0 otherwise É 


The integration between quantum computing and BPSO is known as QBPSO [28]-[30]. The use of 
the contraction-expansion coefficient is one of the parameter differentiators in QBPSO to control the 
convergence velocity from the particles. The first search, which is more global and dynamic, is 
accommodated by using an initial Bmax value of 1. Then, the beta (£) value is gradually decreased until it 
reaches min = 0.4. This is done to complete the QBPSO algorithm search with a better local search. In (18) 
gives the evolution of B. 


BCE) = Bmax — (=m). iter(t) (18) 


itermax 


The new positions in quantum particle swarm optimization (QPSO) are provided by (19) and (20) 
respectively, using the Monte Carlo method. It is noteworthy that time units are used for particles to move. 
Additionally, in the migration process, there is an evaluation process, where the best position is represented 
by mbest. From the QPSO result, the velocity value will be converted to binary by using (16) and (17). 


Xia (t + 1) = Pia (t) + P (t). (mbest g(t) — Xia (t)).ln In ©, Pia (t) 2 0.5 (19) 
Xia (t + 1) = Pia (t) — P (t). (mbest g(t) — Xia (t)).ln In ©, Pia (t) < 0.5 (20) 
With Pia (t) and pa(t) values like (21) and (22). 


Pia (t) = Palt). pbestia (t) + (1 — pa (t)). gbestia (t) (21) 


palt) = tu (22) 


(c1-ria(t)) + (C2-r2a(t)) 


3. RESULTS AND DISCUSSION 

In this research, a 20 kV 34-bus radial network is utilized with the OLTC’s location on the IEEE 34 
bus system depicted in Figure 3. The net and load data on each bus are presented in Tables 1 and 2. The bus 
data table includes a load-type column where the value 1 denotes constant power, 2 denotes constant current, 
and 3 represents constant impedance. In the system, there are 34 buses and five transformers, positioned 
between specific buses. The power flow in the radial distribution system is calculated using the BIBC-BCBV 
method with average lines and loads taken from (13). The total active power of the loads is 428.638399 kW. 
Without OLTC optimization, the power loss from the distribution system is 21.756107 kW, and the average 
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voltage magnitude is 19.00467941 kV, with the lowest voltage at 18.5581 kV. These values are obtained with 
tap formation for transformers | to 5 in a sequence of -1, 0, 2, -1, 0. 


YY 


Duol 


20 kV 
Buo2 
20 kV 
Buo3 
20 kV 
Qn 
Buod CO S WA 
20 kV v 
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20 kv Duos 
nV 
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20 kV 
Buo9 
20 kV 
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Duels Buold 
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20 kv 
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Figure 3. Single line diagram distribution system 34-bus 


The voltage drops and current values in every bus and line obtained from the power flow simulation 
in Table 1 represent the condition or characteristic of the IEEE 34 bus system and are used to calculate the 
power loss. Figure 4 shows the bus voltage of the system within the defined tolerance range of + 5% of the 
base voltage of 20 kV, where buses 16 to 34 experience under voltage. Tap formation for transformer 1 to 
transformer 5 is performed in a sequence of -1, 0, 2, -1, 0. The current at the line between bus 3 and 4 for 
transformer 1 is 20.7 A with tap transformer -1. Similarly, for transformer 2, the current at line 9 between bus 
9 and 13 is 20.7 A with tap transformer 0. For transformer 3, the current at the line between bus 15 and bus 
16 is 20.7 A with tap transformer 2. The current at the line between bus 21 and bus 25 for transformer 4 is 
11.27 A with tap transformer -1, while for transformer 5, the obtained current is 0.57 A with tap 
transformer 0. 


TELKOMNIKA Telecommun Comput El Control, Vol. 22, No. 2, April 2024: 488-498 


TELKOMNIKA Telecommun Comput El Control Oo 493 


Table 1. Power flow result without OLTC optimization using QBPSO 


From bus To bus P ie oa AR) IA) wens (deg) Voltage drops (kV) 
1 2 0.26 0.26 20.07 -33.29 0.02 
2 3 0.18 0.18 20.07 -33.29 0.01 
3 4 3.29 3.28 20.07 -33.29 0.23 
4 5 0 0 0 0 0 
4 6 3.82 3.82 20.07 -33.29 0.26 
6 7 3.03 3.03 20.07 -33.29 0.21 
7 8 0 0 20.07 -33.29 0 
8 9 0.05 0.03 20.07 -33.29 0 
9 10 0 0 0 0 0 
10 11 0 0 0 0 0 
11 12 0 0 0 0 0 
9 13 1.5 1.1 20.07 -33.29 0.09 
13 14 0 0 0 0 0 
13 15 0.65 0.34 20.07 -33.29 0.04 
15 16 3.01 2.24 20.07 -33.29 0.19 
16 17 0.07 0.05 19.55 -33.46 0 
17 18 0 0 0 0 0 
17 19 5.14 3.76 19.55 -33.46 0.33 
19 20 0 0 19.55 -33.46 0 

20 21 0.23 0.17 11.27 -38.58 0.02 
21 22 0 0 0 0 0 
20 23 0 0 8.39 -26.57 0 
23 24 0.19 0.19 8.39 -26.57 0.03 
21 25 0.27 0.2 11.27 -38.58 0.03 
25 26 0.01 0.01 9.32 -38.54 0 
26 27 0.04 0.03 9.32 -38.54 0.01 
27 28 0 0 1.38 -39.2 (0) 
28 29 0 0 1.38 -39.2 0 
25 30 0 0 1.95 -38.82 0 
30 31 0 0 0.57 -37.87 0 
31 32 0 0 0.57 -37.87 0 
31 33 0 0 0 0 0 
33 34 0 0 0 0 0 
Table 2. OLTC optimization result using QBPSO 
Line location Tap transformer 
Initial bus Final bus 
OLTC1 3 4 - 
OLTC2 9 13 
OLTC3 15 16 3 
OLTC4 21 25 -3 
OLTC5 30 31 3 
20.5 
20 
a 195 
2 
= 19 
w 
00 
fa- 18.5 
o 
> 18 
17.5 


1 3 5 7 9 11131517 19 21 23 25 27 29 31 33 


Bus 


Figure 4. Voltage profile without OLTC optimization 


The power loss simulation was performed to determine the active and reactive power losses in the 
transformer lines. The results showed that the active power losses in transformers 1, 2, and 3 were relatively 
high at 3.29 kW, 1.5 kW, and 3.01 kW, respectively, while the losses in transformers 4 and 5 were lower at 
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0.24 kW and 0 kW, respectively. The reactive power losses in all transformers were also recorded, with the 
highest value of 3.28 KVAR in transformer 1 and the lowest value of 0 KVAR in transformer 5. These losses 
were attributed to the large current flowing through the transformers. Furthermore, the voltage drop analysis 
revealed that the lines between bus 3 and bus 4, bus 15 and bus 16, and bus 21 and bus 25 had relatively high 
voltage drops of 0.23 kW, 0.19 kV, and 0.03 kV, respectively. The findings from this study provide valuable 
insights into the power loss and voltage drop characteristics of the IEEE 34 bus system, which could be 
useful for improving the system’s overall efficiency and performance. 

The QBPSO optimization simulation was utilized to obtain optimal OLTC tap coordination to 
minimize power loss and maintain voltage within a predetermined standard. The result of the OLTC 
simulation using QBPSO is presented in Table 2, which brings changes to the power flow in the system as 
shown in Table 3. The power loss obtained from this optimization is 19.157321 kW, with an average voltage 
of 19.93068 kV. Figure 5 indicates that there is no under-voltage bus, and the minimum voltage obtained is 
19.099 kV. The required time for the optimization is 3.858783 seconds, as presented in Table 4. 


Table 3. Power flow result with OLTC optimization using QBPSO 


From bus To bus P ie oy AR) I(A) wee? (deg) Voltage drops (kV) 
1 2 0.25 0.25 19.47 -33.53 0.02 
2 3 0.17 0.17 19.47 -33.53 0.01 
3 4 3.09 3.09 19.47 -33.53 0.22 
4 5 0 0 0 0 0 
4 6 3.6 3.59 19.47 -33.53 0.26 
6 7 2.85 2.85 19.47 -33.53 0.2 
7 8 0 0 19.47 -33.53 0 
8 9 0.04 0.03 19.47 -33.53 0 
9 10 0 0 0 0 0 
10 11 0 0 0 0 0 
11 12 0 0 0 0 0 
9 13 1.14 0.77 19.47 -33.53 0.09 
13 14 0 0 0 0 0 
13 15 0.61 0.32 19.47 -33.53 0.04 
15 16 2.36 1.55 19.47 -33.53 0.18 
16 17 0.07 0.05 18.89 -33.73 0 
17 18 0 0 0 0 0 
17 19 4.8 3.52 18.89 -33.73 0.31 
19 20 0 O 18.89 -33.73 0 

20 21 0.23 0.17 11.37 -38.53 0.03 
21 22 0 0 0 0 0 
20 23 0 0 7.62 -26.57 0 
23 24 0.16 0.16 7.62 -26.57 0.03 
21 25 -0.26 -0.03 11.37 -38.53 0.03 
25 26 0.01 0.01 9.51 -38.48 0 
26 27 0.04 0.03 9.51 -38.48 0.01 
27 28 0 0 1.34 -39.15 0 
28 29 0 0 1.34 -39.15 0 
25 30 0 0 1.86 -38.8 0 
30 31 -0.02 -0.02 0.52 -37.87 0 
31 32 0 0 0.52 -37.87 0 
31 33 0 0 0 0 0 
33 34 0 0 0 0 0 


Voltage (kV) 
© 
uw 


1 3 5 7 9 1113 1517 19 21 23 25 27 29 31 33 
Bus 


Figure 5. Voltage profile with OLTC optimization using QBPSO 


TELKOMNIKA Telecommun Comput El Control, Vol. 22, No. 2, April 2024: 488-498 


TELKOMNIKA Telecommun Comput El Control Oo 495 


Table 4. Comparison of time running program results 


Method Time 
BPSO 10.5559 
QDE 7.182868 

QBPSO 3.858783 


The tap formation for transformers 1 to 5 in sequence is used: -1, 1, 3, -3, 3. The current in the line 
between bus 3 and bus 4 is lowered from 20.7 A to 19.47 A with tap transformer -1 at transformer 1. The 
current in the line between bus 9 and bus 13 is lowered from 20.7 A to 19.47 A with tap transformer 1 at 
transformer 2. The current in the line between bus 15 and bus 16 is lowered from 20.7 A to 19.47 A with tap 
transformer 3 at transformer 3. The current in the line between bus 21 and 25 is increased from 11.27 A to 
11.37 A with tap transformer -3 at transformer 4. The current from 0.57 A to 0.52 A is obtained with tap 
transformer 3 at transformer 35. From the obtained current with the tap transformer, it can be seen that the 
current flowing can be reduced, resulting in a decrease in power losses. 

The simulated power loss in the lines with transformers shows that the active power loss in the line of 
the transformer | decreased from 3.29 kW to 3.09 kW and the reactive power loss decreased from 3.28 kKVAR 
to 3.09 KVAR. The active power loss in the line of transformer 2 decreased from 1.5 kW to 1.14 kW and the 
reactive power loss decreased from 1.1 KVAR to 0.07 KVAR. The active power loss in the line of transformer 
3 decreased from 3.01 kW to 2.36 kW and the reactive power loss decreased from 2.24 KVAR to 1.55 KVAR. 
The active power loss in the line of the transformer 4 decreased from 0.27 kW to -0.26 kW and the reactive 
power loss decreased from 0.2 KVAR to 0.03 kVAR. The active power loss in the line of transformer 5 
decreased from 0 kW to -0.02 kW and the reactive power loss decreased from 0 KVAR to -0.02 KVAR. The 
voltage drops across the lines between bus 3 and bus 4 decreased from 0.23 kV to 0.22 kV, the voltage drops 
in the line between bus 9 and bus 13 remained at 0.09 kV, the voltage drops in the line between bus 15 and 
bus 16 decreased from 0.19 kV to 0.18 kV, the voltage drops in the line between bus 21 and bus 25 remained 
at 0.03 kV, and the voltage drop in the line between bus 30 and bus 31 remained at 0.0 kV. Table 5 presents a 
comparison of simulation results using various artificial intelligence methods. In this paper, GA, BPSO, 
QDE, and QBPSO are compared to QBPSO. Based on Table 5, it is shown that QBPSO converges faster than 
the other methods, as demonstrated in Table 4. In the first experiment, QBPSO was able to directly find the 
minimum power loss of 19.1573 kW, while GA produced a simulation result of 19.3016 kW. 


Table 5. Comparison of running program results 


No GA BPSO QDE QBPSO 
1 19.8232 19.3411 19.3718 19.1573 

2 19.8232 19.3718 19.1573 19.6779 

3 19.8232 19.1573 19.1573 19.4044 

4 19.7020 19.1878 19.1878 19.3718 

5 19.7020 19.2201 19.1573 19.7060 

6 19.3016 19.6484 19.3411 19.4389 

7 19.5169 19.2201 19.1573 19.4044 

8 19.6834 19.1573 19.3411 19.2545 

9 19.7602 19.2201 19.3411 19.4389 
10 19.6828 19.3411 19.1573 19.4279 
Mean 19.6819 19.2865 19.2370 19.4282 
STD 0.162869 0.149271 0.097111 0.165871 
Worst 19.8232 19.6484 19.3718 19.7060 
Best 19,3016 19.1573 19.1573 19.1573 


4. CONCLUSION 

The QBPSO method was employed in this paper to determine the optimal arrangement of tap 
transformers in a network for the purpose of minimizing power loss and improving the average voltage level. 
The tested IEEE 34 bus system was found to have a power loss of 21.756107 kW and an average voltage 
magnitude of 19.00467941 kV, with the lowest voltage recorded at 18.5581 kV. The OLTC values for 
transformers 1 to 5 in sequence were -1, 0, 2, -1, and 0. Through OLTC optimization, QBPSO simulation 
resulted in a power loss of 19.157321 kW and an average voltage of 19.93068 kV, with OLTC values for 
transformers 1 to 5 in the sequence being -1, 1, 3, -3, and 3. 
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